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Abstract 

Nonstationary properties of ballistic constrictions NcS'S with disordered S'S 
electrodes are analyzed theoretically. Amplitudes of Andreev and normal re- 
flections at the constriction are related to the solutions of a stationary Green 
function problem in an inhomogeneous S'S electrode in a dirty limit. This 
provides a generalization of the model of Blonder, Tinkham and Klapwijk for 
a spatially inhomogeneous case. The relation between quasiparticle current in 
NcS'S junctions and energy spectrum of a S'S proximity sandwich is found 
for arbitrary parameters of S' and S materials and of S'S interface. 

PACS numbers: 74.50.+r, 74.80.-g, 74.80.Fp 
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I. INTRODUCTION 



Tunnel junctions with high critical current density are presently a subject of extensive 
experimental investigation (see [|l| and references therein). For such junctions the trans- 
parency of a tunnel barrier is not small, and therefore they do not fulfill the conditions of a 
standard tunnel theory. As was shown in |T[], barriers in Nh — AlOx — Nh junctions with high 
Jc are likely to be a series of constrictions, each having rather large transparency. Therefore, 
a small ballistic ScS constriction is a suitable starting point to discuss more complicated 
models for high Jjunctions. 

The properties of NcS constrictions are well understood in the framework of the model 
of Blonder, Tinkham and Klapwijk (BTK) 0. In their approach, current through a con- 
striction is fully determined by the amplitudes of normal and Andreev reflections at the A^^" 
interface. The model assumes that both N and S metals are in thermal equilibrium and 
that both are spatially homogeneous. Later, the model was further developed by Klapwijk 
et al, and Octavio et al (KBT and OTBK models) to treat ScS constriction, in which 
the subharmonic structure on I — V curves was explained to be due to multiple Andreev 
reflections at the constriction. Whereas the condition of thermal equilibrium is generally 
fulfilled for the constriction geometry, the other condition of spatial homogeneity of the su- 
perconducting electrodes is less general. An important case of an inhomogeneous system is 
the SS'cS'S junction, where S' is a superconductor with < T, (T^ = corresponds tothe 
particular case of a normal metal). 

Previous work on the generalization of the BTK approach to account for spatial inhomo- 
geneity was started by Van Son et al P]. Andreev reflection in the NcN'S system (where N' 
is a normal metal) was considered for a gradual variation of a pair potential near the normal 
metal - superconductor interface. To model Andreev reflection from the N' region with 
proximity induced superconductivity, the existence of a spatially dependent pair potential 
A„/(x) was assumed in A^'. However, whereas the Cooper pair density in a normal region 
is indeed nonzero due to the proximity effect 0, the pair potential A„/(a;) = 0, if = 0. 
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Therefore a more consistent approach is needed to modify the Andreev amphtudes in the 
NcN'S sandwich in comparison to the BTK case of NcS. Such a theory is necessary, in par- 
ticularly, to interpret experiments on point contact spectroscopy of proximity systems, like 
the one performed recently on a bilayer consisting of doped Si backed with superconducting 

iV6[g. 

It will be shown in this paper that the existence of a pair potential A„/(x) in A^' is 
not a necessary condition for Andreev reflection at the NcN' boundary. As is known, the 
spectrum in a normal region is modified due to the proximity effect, namely, bound states 
exist in A^' layer at energies below the energy gap of a superconductor, [1,0. The energy 
of the lowest bound state corresponds to an energy gap A^^/ induced in A^', which for a thin 
enough layer is close to A^. Therefore it is clear qualitatively that in NcN'S systems Andreev 
reflection processes take place not only at the N'S boundary, but also at the NcN' boundary, 
because at E < Agn> a quasiparticle can not penetrate into A^'. This result however is not 
directly evident from Bogolubov de Gennes (BdG) equations, where the pair potential A„/(x) 
plays the role of an effective scattering potential for nondiagonal scattering. To find such 
a potential for N'S sandwiches, one should either solve a complete three-layer problem in 
the BdG equations, taking into account both A^A^' and N'S boundaries simultaneously, 
or to use the microscopic Green functions approach. For a disordered S'S system the first 
approach would require knowledge of the full scattering matrix, while the second one is much 
more straightforward, and we shall use it in this paper to find the nondiagonal potential 
for tunneling into the S'S bilayer and to calculate the quasiparticle current for the NcS'S 
contact. 

Previously, a microscopic approach based on Green functions formalism was used to 
study the properties of NcS and NcN'S microcontacts without impurity scattering (clean 



limit). For the NcS case, Zaitsev |]10[ has derived boundary conditions for the quasiclassical 
Eilenberger equations at the contact interface and has calculated the current through the 
contact at arbitrary transparency of the interface, thus providing a microscopic derivation 
of BTK model. Independently, the same derivation was done by Arnold with another 



method of retarded Green functions. As a particular case, properties of clean NIN'S con- 
tacts with arbitrary barrier transparency were considered in |Tl|]. More recently, detailed 
calculations for the NcN'S case were done in the framework of the Arnold model in 
with the additional assumption of small transparency barrier between A^' and S layers. 

In a case of strong disorder (dirty limit) nonequilibrium aspects got a lot of attention 
recently. In [[ISHT^ the conductivity of a dirty NN'S contact was studied in the model, where 



A^ and S electrodes are reservuars with fixed electrical potentials whereas the disordered 
contact region A^' is a one-dimensional bridge in nonequilibrium. Such a situation can be 
realized in disordered NS contacts when the A^' layer simulates the interface region of the 
order of the inelastic mean free path. An enhancement of the zero bias conductance a{0) 
was predicted in this system. 

In the present paper we consider a physically different situation: a ballistic constriction 
NcS'S of the size smaller than the mean free path having disordered electrodes. In this case 
the potential drop takes place at the constriction, and therefore the electrodes are in thermal 
equilibrium. The ballistic condition is important. It should be compared with the opposite 
limit of disordered constriction of a size larger than a mean free path. The latter case was 
studied theoretically by Artemenko, Volkov and Zaitsev for ScS constrictions and more 



recently by Volkov [|I^ for NcN'S constrictions with large transparency. It was obtained, 
in particularly, that due to disorder no conductance doubling is present at zero bias V = 0. 
As is shown in this paper, the conductance of the ballistic NcN'S constriction is different 
in many respects as zero bias conductance doubling is being present for large transparency. 
The influence of the finite transparency of a constriction is also studied. Physically, our 
approach is closely related to the BTK one. We will give expressions for coefficients of 
Andreev and ordinary electron reflections at the ballistic constriction and their relation to 
the energy spectrum of the disordered SS' system that is investigated in terms of Green 
functions. As a result, the simplicity and direct physical meaning of the BTK solutions are 
combined with the approach based on the selfconsistent solution of a stationary dirty limit 
Green functions problem. 



II. THE MODEL 



Let us consider the boundary between N and S'S as a small constriction of size a <S 
min(/„, Ig'), where Ig' are mean free paths of N and S'. The constriction is characterized 
by a transmission coefficient 



D 



(1) 



where H is strength of repulsive potential Hd{x) located at the NS' interface and vix,V2x 
are the components of the Fermi velocities of N and S' normal to the interface, respectively. 
We assume that S' and S metals are in the dirty limit /^/^ <^ 3, whereas in A'^ there is no 
limitation on the mean free path. 

Let us follow the BTK notations. In the BdG equation formalism, the excitations are 



represented by a vector ip 



9[x) 



, where f{x) describes electron- like excitations in a 



superconductor, andg{x) describes hole-like (time-reversed) excitations. The BdG equations 
have the form: 



"''dt 



2m 



- fi{x) + V{x)) fix) + A{x)g{x, t) 



(2) 



^^S = - (-Sf - /^(^) + n^)) 9{x) + A(x)/(x, t) 
where ii{x),/\{x) and V{x) are the electrochemical potential, the pair potential, and the 
ordinary potential, respectively. 

In the absence of impurity scattering, the transmission and reflection coefficients for a 
quasiparticle incident from clean N to clean S, are found in the BTK model by considering 
incoming, reflected and transmitted waves near the NS boundary: 

(A 



(3) 
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2--(l + Ve2-AVe) 



(5) 



where A is a the bulk energy gap of S", mi 2 are effective masses of the contacting metals 
and e is the quasiparticle energy. 

Matching these solutions at the NS boundary, one can find the energy- dependent Andreev 
reflection coefficient A{e) and normal reflection coefficient B[e): 

A{e) = aa* - 

B{e) = bb* 



- (l+Z2(l-|,,fjj2' 



(6) 



(l+(l-|r7|^)Z2)2 



^0 e 



(7) 



Uo{e) 1 + E/Ve^ - 

where Z is related to normal transmission coefficient D by: (1 + Z"^) = D^^. It is important 
to note that both A{e) and B{e) are controlled by one energy- dependent parameter 77(e), 

I 1 2 I 1 2 

\ri\ being simply the ratio of probabilities for the excitation to be in hole-like, |fo| , or in 
electron-like, ImqI^, states. 

The solutions (5) for Uo,Vo are applicable for a clean homogeneous BCS superconductor, 
and they are not valid in spatially inhomogeneous dirty S'S sandwiches. To find the coeffi- 
cients A{e) and B{e) in the latter case one needs to calculate 77(e) as a solution of Gor'kov 
equations (GE) in the SS' system. For the ballistically clean constriction we can take the 
advantage that at a length scale smaller than the electron mean free path the solution of 
the GE can be written in the form of a combination of plane waves. 

The GE in the region near the constriction, S', have the form [|T^: 

{2m2(e'(a;) + l^) + £2} G,{x, x') + 2m2A'{x)F:{x, x') = 6{x - x') 
{27n2(e'(x) - ;x) - ^} F,(x, x') + 2m2A'(x)G,(x, x') = 



Here E' and A'(x) are the energy and pair potential renormalized by impurity scattering 
according to 



e\x) = e + i {G,{x, x)) /2r, A'(x) = A(a;) + x)) /2r, 



(9) 



T = (27rcy^A^(0))~^ is scattering time, c and V are impurity concentration and scattering 
potential, respectively. The brackets (...) denote angle averaging, and Ge{x,x') and F^{x,x') 
are the normal and anomalous Green functions in energy representation. 
The pair potential in (8) and (9) is given by the selfconsistency equation: 



A{x) = gTj2Fe=-iw„{x,x) 



(10) 



where g is the coupling constant and Un = 7rT(2n + 1) is the Matsubara frequency. Note 
that in the particular case of = the pair potential in S' is zero, A{x) = 0, whereas 
F^{x, x) is finite. 

We note that {G^{x, x)) and {F^{x, x)) are quasiclassical Green functions of a dirty super- 
conductor. They obey diffusion-like equations PUipTI with the boundary conditions derived 



m 



Let us consider the solutions of the GE in the S' region at distances less than C,s' from 
the constriction. In this region one can neglect variations of A'{x) and e'{x) and to write 
the solutions of the linearized equations as a combination of plane waves: 
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= C{x') 













-iq2 X 



(11) 



where the slowly varying functions g{x), f{x) determine amplitudes of electron like and hole 
like excitations. Substitution of the solution eq.(ll) to the GE leads to the following linear 
system of equations for g{x) and f{x): 



^ 9{x) ^ ^ 



\ fix) J 



2m2 [e + i{G,{x)) /2t + ^i]- ql [A + (F, (x) ) /2r] 

2m2 [A + /2r] [e + i /2r - /i] + ql j 



^ ^ 9{x) ^ 



fix) ^ 
(12) 



Here is determined by the dispersion relation DetL — 0, which in the dirty hmit At <^ 1 
leads to the following result 



2m2 



2r 



(13) 



Eq.(ll) describes transmitted electron and transmitted hole waves in the S' region near 
the constriction. Note that, in general, reflected electron and reflected hole waves are also 
present in S' due to impurity scattering. However, the condition of a ballistic constriction, 
a <^ Is', made it possible to neglect these waves in eq.(ll), because under this condition 
waves in S' are scattered diffusively away from the contact, and there is a small probability 
for a wave, scattered at the distance Ig' from the constriction, to reach it again. 

To find transmission and refiection probabilities for a given problem, one should compare 
the solutions (4) and (11) at distances smaller than mean free path Igi from the constriction. 
In this region one can neglect the small energy terms in (3) and (13) in comparison with 



chemical potential /i, i.e. set q'2 = = v^^m^ the phases of all transmitted waves. 
Then, to find the probabilities A(e) and B{e) one should find the ratio f{x)/g{x). For the 
dirty limit we obtain from eq.(12): 

g{x) l + {G,{x))- ^ > 

Finally, substituting 77(e) in eqs.(6),(7) by f{x)/g{x) from eq.(14) and using the normal- 
ization condition for the Green functions {Gf) + {F'^) = 1, we find that the Andreev and 
normal refiection coefficients are directly related to the local energy spectrum of S' near the 
constriction: 

m = „ (15) 

^' |l + 2Z2+(G,(0+))|" ^ ' 

""^'^ ^ |l + 2Z^ + (G,(0+))r ^''^ 
where (Gg(0-|-)) and (^^(0+)) are quasiclassical Green functions in the vicinity of the contact 
in S' . The expressions (15) and (16) are the central result of the paper. They generalize the 



corresponding BTK relations p| for a spatially inhomogeneous case. For the latter case the 
Green functions are given by {G^{0+)) = —ie/ Aq — e^, {F^{0+)) = Aq/^J Aq — e"^ , and it 
is easy to check, that the BTK relations follow from eqs.(15), (16). 

The coefficients A{e) and B{e) given by eqs.(15,16) fully determine the quasiparticle 
current through NcS'S contact. Namely under the condition of thermal equilibrium in both 
electrodes the current through the constriction is expressed via A{e) and B{e) in complete 
analogy with the BTK result for the NcS constriction |0| : 

m = -^^f^ [fo{e + eV) - /o(e)] [1 + A{e) - B{e)] de (17) 

where Rq = [2Ni{0)Se'^vpi] ^ is the Sharvin resistance, 5* is contact area, A''i(0) and vpi 
are the density of states per spin and Fermi velocity of electrode, respectively (we have 
taken vpi < vf2)- fo is the Fermi distribution function and the brackets (...) denote angular 
averaging. The relations (15) and (16) show, that in the dirty limit local information is 
given by point contact measurements. In particular, the local density of states near the 
constriction is given in the usual way as A^(e) = Re {{Ge{0+))}- It is seen however from 
eqs.(15)-(17), that generally (for arbitrary Z) the current is not determined solely by A^(e) 
like in a tunnel theory, but rather a crossover to the limit of a tunnel NIS'S junction takes 
place for small transparencies, like in NcS contacts 0. We note, that the case of SS'cS'S 
ballistic constrictions can be considered in the same way, by extension of the OTBK model 
§ with the help of eqs.(15), (16) for A{e) and 5(e) |2|. 

Therefore, the problem is reduced to calculation of the functions {Ge{0+)) and (Fe(0+)) 
for the dirty SS' sandwich. 

III. PROXIMITY EFFECT IN THE DIRTY SS' SANDWICH 

The proximity effect in dirty S'S sandwiches was studied previously in [0,0 for the 
case of arbitrary transparency of the S'/S interface, where the case of a thin 5" layer was 
considered. Now we first generalize the results of p2| , P^ to arbitrary S' layer thickness and 
then, using these solutions, calculate the quasiparticle current of a NcS'S contact. 
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The angle averaged quasiclassical Green functions in the dirty S'S bilayer satisfy the 
equation: 

ieFs>^s{x) + \Gs',s{x) — Fs'^s{x) j - As> ,s{x)G s' ,s{x) = 0, (18) 

where we have omitted the angular brackets (...) for the functions Gs'.s, -^s',s in S' and S 
regions respectively. Here Dg/ g is diffusion coefficient, A^'^^ is the order parameter. To be 
more specific, we shall discuss below the particular case of A^' = 0, i.e. of Tc' = (A^^" 
sandwich). The generalization to the case of nonzero Tc' is straightforward ||2^ and does 
not change our results qualitatively. The latter case will be discussed separately elsewhere. 

It is convenient to rewrite eq.(18) using the notations G{e,x) = cos6{e,x), F{e,x) = 
sm6{e,x). Then eq.(18) takes the following form in and 5* regions: 

CN,s^N,si^) + sin 9N,six) + An,s{x) cos 6N,six) = 0, (19) 
with the boundary conditions at the NS interface (x = 0) ^ 



in the bulk of the S"— layer 



A,(x)ln^ + 2^^ 



A.(x) 

sm tls[x, e = tuin) 



The parameters ■jb and 7 



(20) 



9s{oo) = arctan(i Ao (T) /e) , (21) 
as well as at the metal free surface {x = —(1^): 

e'^i-dN) = 0. (22) 
The selfconsistency equation for the order parameter in S region has the form: 

: 0. (23) 



Rb PsCs /^.x 

Pn<;n Pn<;n 
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have simple physical meanings: 7 is a measure of the strength of the proximity effect be- 
tween the S and metals, whereas 7^ describes the effect of the boundary transparency 



between these layers. Here pN,s, ^n,s = y Dn,s/'^'^Tc and Dn,s are normal state resistivities, 
coherence lengths and diffusion constants of and S metals, respectively, while Rb is the 
product of the resistance of the NS boundary and its area. We have normalized e and A{x) 
to vrTc , where Tc is the critical temperature of the bulk S. 

Previously selfconsistent solutions of the boundary value problem (19)- (24) for the NS 
sandwich at arbitrary values of 7 and 7^ were studied only for a thin layer, dj^ ^ C,n, in 
2^^6|. In particular, the densities of states in N and 5* layers N^^si^) = -Re (cos6'iv,s(e)) 



were discussed. It was shown for this case that for any values of 7 and 75 a superconducting 
state is induced with finite energy gap Ag^. When both layers are thin dN,s ^n,s, the 



results of the well known McMillan tunnel model |27] of the proximity effect can be repro- 



duced |2^. Then, in the language of the McMillan model, the existence of a finite gap in the 



considered diffusive limit is due to finite inverse lifetime of quasiparticles in thin layer. 



N 



.-1 



{hvpN /2dN)D, where D is the angle averaged NS boundary transparency. 



As a result, there is a finite minimal energy of quasibound states in A^ and therefore, a 
nonzero energy gap is induced. For finite thickness d^ the densities of states in a dirty NS 
sandwich were discussed previously only in the framework of rigid boundary conditions, i.e. 



for 7/75 < 1, in IIMl- 

In the general case of arbitrary dj\f,ds,'y and •jb the solutions should be determined 
selfconsistently together with the spatial dependence of the order parameter from the self- 
consistency equation. The results depend essentially on the parameters 7, 7b, but from 
lifetime considerations it is clear qualitatively, that the above conclusion about a finite gap 
should hold also for any finite thickness of A^. Below we will calculate the density of states 
in A^ from the solution of eqs.(19)-(24) and will demonstrate the reduction of the gap in A^ 
with increase of dN- 

Taking advantage of the condition A^r = one can integrate the eq.(19) in A^ region and 
with the help of boundary conditions (20) obtain 
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COS^Ar(O) — COS^Ar(— dAf) 



sin^(gg(0)-g^(0)) 



(25) 



The analytic solution of eq.(25) is simplified only in a limiting case of small 7/7^ ratio, 
when one can take 9s{0) — ^s(oo) in the first approximation. Then for sufficiently low 



we obtain the BCS expression with a gap Aq for the density of states in this energy range: 



7V]v(e, 0) = Re ye/ye^ — AqJ . For large N layer thickness, ^jv, one can substitute 

cos9N{d) = 1 in eq.(25) and find the asymptotic behavior of A^Ar(e) at the NS boundary: 



For the case of arbitrary N layer thickness and arbitrary values of the parameters 7, 7^ 
the boundary value problem (19)-(24) was solved numerically. The results of the calculations 
of the densities of states in at x = and x = — djv at low temperatures T <C Tc are 
presented in Figs.l and 2 for a number of d^/^N ratios. As is seen from Fig.l aX, x — —dN 
(free surface of N) two peaks exist in A^(e) for small 7 values, provided that dj^/^N < 1. The 
first peak is at e = A^jv and the second one at e = Aq. It is also seen from comparison of 
Figs.l and 2, that the second peak at e = Aq is smeared out with the increase of c^at as well 
as with the increase of 7. On the other hand, the peak at e = A^tv becomes more pronounced 
at large djsi- Dotted lines show the behavior of N{e) in N at the boundary with S {x = Q). 
The asymptotic behavior given by eq.(26) should take place at large d^- It is important 
to note, that the energy gap A^^r is preserved for all djq, going to zero rather slowly. This 
is consistent with the qualitative picture of the gap being proportional to inverse lifetime 
in N ^ which in the considered diffusive approximation is given by r^^ ~ {hDN/d%)D. The 
results of study of the dependence of A^jv on the parameters of the NS sandwich will be 
presented in more detail elsewhere. 

Using the solutions 9n{€, —d^), one can calculate the reflection coefficients A{e) and 
S(e) and then the quasiparticle current for the NcNS contact from eqs.(15)- (17). To 
be more specific, let us discuss here the case of a thin N layer, d^/^N ^ 1- Then the 



energies e <^ 7^^, it follows from (25) that ^Ar(O) ~ ^5(0) = arctan(iAo/e), and as a result 





(26) 



12 



number of parameters is reduced from 7,75 and to the following set: 7^ = •ydN/^N and 



7bn = IsdN/iN 1231 • 

Fig. 3 shows the results of the calculations of the reflection coefficients A{e) and B{e) for 
a number of Z values. It is seen that for Z 7^ a characteristic two-peak structure exists for 
both A(e) and B{e), which is directly related to the two-peak structure of Njq{e) discussed 
above, the first peak being at e = A^jv and the second one at e = Aq. 

The zero-temperature conductance of NcNS junctions calculated according to eq.(17) 
is shown in Fig.4 for the same parameters. Again the two-peak structure is present in 
dI{V)/dV at voltages AgN and Aq. In accordance with the arguments given above, the 
peak at Aq would be smeared out quite easily by any pair-breaking process, i.e. by large 
7m values (spatial gradients in S), or large d^- Then such a spectroscopy will show only 
a proximity induced energy gap in with almost no signatures of Aq. The position of the 
first peak at eV = Ag^ can be used to study properties of the NS interface for any given 
material combination. The appearance of such a conductance peak at low bias for Z 7^ is 
a consequence of the given model, which follows directly from the structure of the densities 
of states in the region at x = —dj^j, as is seen from Figs.l and 2. 



IV. CONCLUSIONS 

In conclusion, the BTK model is generalized for a spatially inhomogeneous case of NcS'S 
ballistic constrictions with disordered S'S electrodes. The expressions for the amplitudes 
of Andreev and normal refiection are given, which allow to calculate a quasiparticle current 
for arbitrary parameters of S' and S materials and their interface, if the conditions of the 
dirty limit are fulfilled. An energy gap in S' is always present, even for finite thickness of 
the S' layer. The magnitude of this gap is studied as a function of the parameters of the S 
and S' materials, as well as of the transparency of the SS' interface. It is shown, that the 
conductance of ballistic NcS'S junctions refiects proximity induced energy gap in S' and, 
under certain conditions, also the bulk gap of the superconductor S. 
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FIGURES 

FIG. 1. The densities of states in the N' layer of the N'S sandwich, normahzed to their 

normal-state values, at the free surface (solid lines) and at the N'S boundary (dashed lines) for 
different thicknesses (In/Cn = 10 (curve 1), 2 (2), 1 (3) and 0.5 (4). 

FIG. 2. The densities of states in the N' layer of the N'S sandwich, normalized to their 
normal-state values, at the free surface (solid lines) and at the N'S boundary (dashed lines). 
Thicknesses cIn/^n are the same as in Fig.l. 

FIG. 3. Probabilities of Andreev reflection, ^(e) (solid lines), and of normal reflection, B(e) 
(dashed lines), for the ballistic NcN'S constriction with 7^ = 0.1 and jbn = 1- 

FIG. 4. Zero-temperature conductance of the ballistic NcN'S constriction with the same 
parameters as in Fig.3. 
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